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Abstract 

We advocate the use of Daubechies wavelets as a basis for treating a 
variety of problems in quantum field theory. This basis has both natural 
large volume and short distance cutoffs, has natural partitions of unity, 
and the basis functions are all related to the fixed point of a linear renor- 
malization group equation. 



1 Introduction 

Daubechies wavelets and their associated scaling function are functions [JJ [5] [3] 
[3] [5] [6] used in signal processing and data compression. The functions have 
useful properties that have advantages for applications to problems in quantum 
field theory. Daubechies wavelets have been used in some field-theoretic appli- 
cations [2] [8] [9] [10] and different types of wavelets have been used in more 
formal treatments of field theory [T2] [T3] [TTJ [TJ] , however many useful prop- 
erties of wavelets do not appear to have been exploited. Some of the properties 
of Daubechies wavelets that are potentially useful for field theory applications 
are: 

1 . The functions are an orthonormal basis for square integrable functions on 
the line. This means that one can formally express fields exactly as infinite 
linear combinations of discrete field operators and coefficient functions. 
This can be used study the nature of corrections to truncated theories or 
to compute the singular coefficients in operator product expansions. 

2. The basis functions have compact support. This means that discrete fields 
associated with basis functions with space-like separated support com- 
mute. The maximum size of the support of the basis functions can be 
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controlled. There are also infinite numbers of basis functions with arbi- 
trarily small support so it is possible to use these operators study locality 
and how locality is violated in truncations. In these representations local 
fields can be replaced by a local algebra of well-defined operators. 

3. The Daubechies basis functions are related to solutions of a linear renor- 
malization group equation. It is a natural basis to formulate renormaliza- 
tion group transformations. The basis allows for a separation of scales; 
a fine resolution truncated Lagrangian or Hamiltonian can be exactly ex- 
pressed in terms of a coarse resolution truncated Lagrangian or Hamilto- 
nian of the same form, plus some addition operators that fill in the missing 
fine-scale physics. Eliminating the fine scale operators gives a new Hamil- 
tonian, involving the same masses and coupling constants, in the coarse 
scale degrees of freedom that includes the effects of the eliminated fine 
scale degree of freedom. This can be repeated to include the effects of 
eliminated degrees of freedom on arbitrarily small scales. 

4. The basis is generated from a single function, which is the solution of the 
linear renormalization group equation, using discrete unitary translations 
and discrete unitary scale transformations. While the basis functions have 
a fractal structure, overlap integrals of products of any number of the basis 
functions and low-order derivatives can be computed analytically using 
only the renormalization group equation and a scale fixing condition. In 
addition, integrals of products of basis functions and polynomials can also 
be computed exactly using the renormalization group equations. 

5. Subsets of the basis functions contain compact, locally finite partitions of 
unity. This property can be used to investigate symmetries like Poincare 
invariance. The partitions of unity can be used to make local realizations 
of the Lie algebra. The generators can be expressed exactly as sums of lo- 
calized operators obtained by integrating the operator densities multiplied 
by the functions in the partition of unity. 

6. The basis functions have a limited amount of smoothness, and can locally 
pointwise represent low-degree polynomials. This means that overlaps of 
low-order derivatives of field operators can be computed directly, rather 
than by using a finite difference approximation. 

7. Because the basis functions have compact support and there are basis 
functions with arbitrarily small support, the basis has both natural short- 
distance and large- volume cutoffs. 

The purpose of this paper is to discuss some problems in field theory where 
wavelets might have advantages as well as some of the specialized methods for 
dealing with fractal basis functions. In section two we discuss the construction 
of the wavelet basis and introduce the renormalization group equation for the 
Daubechies scaling function. In section three we introduce our free-field conven- 
tions. In section four we introduce the wavelet representation of free fields. This 
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representation is also relevant in interacting theories, when the natural volume 
and resolution truncations in the wavelet representation are used to make finite 
number of degree of freedom truncations of the interacting theory. In section 
five we discuss the scaling properties of the wavelet fields. In section six we dis- 
cuss Poincare invariance and partitions of unity. In section seven we discuss the 
formulation of renormalization group equations in the wavelet basis. In section 
eight we discuss gauge transformations in the wavelet representation. In the 
appendix we compute the overlap coefficients that couple different scales in the 
free-field Hamiltonian. 



2 Basis construction 

The Daubechies scaling function is the solution of a linear renormalization group 
equation subject to a scale-fixing normalization condition. All basis functions 
are generated from the scaling function using discrete translations, dyadic scale 
transformations, and sums. 

The discrete translation operator, T, and a dyadic scale transformation, D 
are unitary operators on L 2 (1R) defined by 

(Df)(x) = V2f(2x) (Tf)(x) = f(x - 1). (1) 

The operator D decreases the support of a function with compact support by a 
factor of two in a manner that preserves the L 2 (R) norm. The scaling function, 
s(x), is the solution of the linear renormalization group equation 

2K-1 

s( x ) = ^ hiDT l s{x). (2) 
;=o 

Equation ^ states that the scaling function is a linear combination of 2K- 
translated copies of itself on a scale smaller by a factor of 2. Since this is a 
homogeneous equation for s(x), we are free to assign the normalization. It is 
fixed by the scale-fixing condition 

J s{x)dx = 1. (3) 

Equations ([2][3|) define the scaling function. Equation ([2]) is called the scaling 
equation. The h n are real coefficients that are characteristic of the type of 
wavelet. The coefficients for some of the Daubechies scaling functions are given 
in Table I. 

Given the solution of ((2][3|), scale-l/2 fc scaling functions are defined 

s k n (x) ~ D k T n s(x) (4) 

by applying n unit translations followed by k dyadic scale transformations to the 
original scaling function. Linear combinations of these functions with squarc- 
summable coefficients span a subspace T-Lk of L 2 (R), which we call the resolution- 
l/2 k subspace. The scaling equation implies the relation 

n k+1 d n k (5) 
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Figure 1: Daubechies K — 3 scaling function. 



and more generally 

Hk+m 3 U k (6) 

for m > 0. This means the lower resolution spaces are linear subspaces of the 
higher resolution spaces. 

The mother wavelet w(x) is the following linear combination of the s^x) 
scaling functions: 

2K-1 2K-1 

w(x) = ]T giDT l s{x) = Yl 9is}{x) (7) 

1=0 1=0 

where the coefficients gi are related to the coefficients hi by reversing the order 
and alternating the signs 

9i = (-) l h 2K -i-i. (8) 

Plots of the Daubechies K = 3 scaling function and mother wavelet are given 
in figures 1 and 2. Scalc-l/2' £ wavelets, w*(x), are constructed by applying 
translations and dyadic scale transformations to the mother wavelet 

w k n {x) = D k T n w{x). (9) 

The scale l/2 k wavelets are designed to span the orthogonal complement, Wfc, 
of H k in ! 

nk+i=n k (BWk. (10) 

The coefficients hi are chosen so the translates of the scaling function and mother 
wavelet are orthonormal and finite linear combinations of translations of the 
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Figure 2: Daubechies K = 3 mother wavelet 



scaling function can be used to make local pointwise representations of polyno- 
mials of degree K — 1. 

Continuing the decomposition (ITU|) inductively leads to the following orthog- 
onal direct sum decomposition of the Hilbert space of square integrable functions 
on the line 



This gives a decomposition of the Hilbert space into orthogonal subspaces with 
successively finer resolutions. Specifically the basis for can be used to repre- 
sent features down to scale l/2 fc , Wk represents features down to scale l/2 fc+1 
that cannot be represented on scale l/2 fc , Wk+i represents features down to 
scale l/2 fc+2 that cannot be represented on scale l/2 fe+1 , etc.. Smearing local 
fields with these basis functions gives discrete fields that are sensitive to average 
information on different scales. 

The coefficients h n , which characterize the Daubechies K- wavelets, are so- 
lutions of the system of equations 



i 2 (R) = U k ®W h ®W k+ i@--- . 



(11) 
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Table 1: 



The first equation is a necessary condition for the scaling equation to have a 
solution; it can be derived by taking the Fourier transform of both sides of 
the scaling equation. The second equation ensures that integer translates of 
the scaling function are orthonormal. The third equation ensures that linear 
combinations of integer translates of the wavelet functions are orthogonal to 
degree K — 1 polynomials. 

This last condition, along with (jlll) . implies that polynomials of degree K—l 
can be locally pointwise represented by locally finite linear combinations of 
scaling functions. 

The solution of these equations for K = 1,2,3 are given in Table 1. The 
solutions are unique up to reversal of order; k — > I = IK — 1 — k. The numerical 
values of the K = 3 coefficients, rounded to two decimal places, are ho = —55 = 
.43, h\ = ,g 4 = 1.11, h 2 = -.93 = -66, h 3 = g 2 = -.34, h 4 = -gi = -.39, h 5 = 
go = —.66. An important observation is that only one of the six coefficients, 
hi = gi has magnitude larger than 1. 

In this work we advocate using the Daubechies K = 3 wavelets. This is 
because of they have small support and have one continuous derivative. The 
derivative allows one to replace the finite difference approximations of deriva- 
tives that appear in the Hamiltonian or Lagrangian with actual derivatives. 

The scaling function s(x) and mother wavelet w(x) have support on the 
interval [0, 2K — 1], It follows from (fTT|) that for any fixed scale, l/2 fc , the 
functions 

{s k n (x)}%L U {<(x)}^oo,™ =fc (15) 

are an orthonormal basis for L 2 (M). The basis functions s„(x) and w^(x) have 
compact support on the interval [2~ k n, 2r k (n + 2K — 1)]. This support justifies 
calling Hk the resolution-l/2 fe subspace. 

The basis functions satisfy the orthonormality conditions 

J dxs k n (x)s k m (x) = 5 mn (16) 
J dxs k n (x)w k m +l (x)= (Z>0) (17) 
J dxw k (x)w l m (x) = S mn S k i- (18) 
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Scaling functions with different k are not orthogonal and are not members of 
the same basis. Similarly, the w^ l (x) with negative I are not orthogonal to 
Sm(x) and are not members of the basis that includes the sj,(i). In general, 
the basis (|15l) consists of scaling functions of a fixed largest scale, 1/2 , and 
wavelets that can represent the structure on all scales l/2 k+l with I > 0. 

A basis for the resolution l/2 fe+m subspace, Wk+m, that includes the basis 
functions for the resolution l/2 fc subspace, Hk, is 

sparx({4+ ro }- _ x ) = spantf**}" U K}^"^)- (19) 

This means the wavelets on scales l/2 fe , l/2 fc+1 , • • • , i/2 fe + m - 1 fill in all of the 
missing higher resolution information in Hk+m that is not contained in T-Lk- 
The scaling functions form a partition of unity 

oo oo oo 

1= J2 s(x-n)= Y, (T n s)(x)= J2 S IW- (2°) 

n=— oo n— — oo n— — oo 

This partition of unity is locally finite, meaning that at any point only a fi- 
nite number of the functions are non-zero, and the functions all have compact 
support. There are also partitions of unity for the scale- l/2 fe scaling functions 

oo ^ oo 

J2 ~m s n^)= E *(2 fc *-n) = l (21) 

n=— oo n=— oo 

except they must be multiplied by the factor -p=. 

The scale- l/2 fe_1 scaling functions and scale- l/2 fc_1 wavelets are both linear 
combinations of the scale l/2 fc scaling functions 



2K-1 

- E h l S 2n + l(*) (22) 
(=0 



2K-1 

w k -\x) = ]T 9is k 2n+l (x). (23) 

1=0 

The inverse of these relations express the scale l/2 k scaling function as a linear 
combination of the scale l/2 fe_1 scaling functions and wavelets 

Sn( X ) = E 71 ™- 2 ™^" 1 ^) +J29n-2mW^ 1 (x). (24) 

m m 

Equations (|22H23[) define the wavelet transform in signal processing. It is a 
pair of filters that decompose a signal into high and low-resolution parts in a 
manner that can be inverted to reconstruct the full signal. In signal procession 
applications it is applied many times to generate multi-resolution decomposition 
of the signal. Equation (l24l) is relevant in field-theory applications, where we see 
that the wavelets represent the lost information when one eliminates small-scale 
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fluctuations. The same coefficients appear in (|22H23I) and (IM1) because these 
two sets of orthonormal basis functions on T-L k are related by a real orthogonal 
transformation. 

For some applications it is useful to express equations (|2"21 and (1231) in infinite 
matrix form 

2n+2K— 1 



where 



(l)= £ (25) 
2n+2K-l 

l (a) = ^ G nm s k m {x) (26) 

m— 2n 

— h"m — 2n G nm — Qrri—2n- (27) 



In this notation becomes 

s n( x ) = ^nrn s m 1 ( x ) + ^ ^nm^m 1 ( X ) (28) 

m m 

where 

^nm — H mn h n — 2m ^~*nm G mn 9n—2m' (29) 

While the matrices are in principle infinite, for fixed n there are only 2K non- 
vanishing terms in the sum, which is clear in the equivalent expressions (|22j) 
and (f23| . 

Note that if we set x = n in the scaling equation it becomes 

2K-X 2n 

s(n) = V2his(2n- 1) = y/2 ^ h 2n - m s{m) (30) 

1=0 m=2n-2K+l 

which has the form of an eigenvalue equation for the matrix M nm — /^n-mi 
where the eigenvectors are (0, s(l), s(2), • • • , s(2if — 2), 0) and the eigenvalue is 
l/\/2- If we differentiate the scaling equation and again evaluate the derivative 
of the scaling functions at integer values we get 

2K-1 2n 

s'(n) = Y 2V2his'(2n- I) = 2^2 ^ h 2n _ m s'{m). (31) 
1=0 m=2n-2K+l 

This is an eigenvalue equation for the same matrix, with eigenvector (s'(0), s'(l), 
s'(2), • • • , s'(2JT - 2), s'(2if - 1)) and eigenvalue 1/2 3 / 2 . Each derivative brings 
down an additional factor of 1/2 in the eigenvalue. Since a finite matrix has 
a finite number of eigenvalues, the scaling function can only be differentiated 
a finite number of times, (n times if the matrix above has an eigenvalues of 
2^(™+ 1 / 2 )). Thus we have the interesting property that the number of derivatives 
is associated with the eigenvalue spectrum of a finite-dimensional matrix. This 
makes it clear that the individual wavelets and scaling functions have only a 
finite number of derivatives. It is interesting that certain linear combinations 
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of scaling functions may have more derivatives than the functions themselves; 
this is the case with the partition of unity (l20l) . The sum is constant, which is 
infinitely differentiable, but the individual functions in the sum are not! 

For fields that are functions of four space-time variables, we use a basis for 
the square integrable functions of the three space variables consisting of prod- 
ucts of the wavelet basis functions, (fTijj) . in each of the three variables. We 
identify basis functions that are the product of three resolution 1/2*-' scaling 
functions as the resolution l/2 fc basis; the seven other products of combinations 
of wavelets and scaling functions span the orthogonal complement of the reso- 
lution l/2 fc subspace in L 2 (R 3 ). We introduce a compact notation to make the 
three-dimensional case look like the one-dimensional case. 

In three dimensions we define: 

4(x) = 4(^)4 (x^ (z 3 ) (32) 

and use the notation to w™ a (x) to represent the remaining basis functions that 
have one of the following seven forms 

<i, fc3 W := 4(xi)4(x 2 X(x 3 ) (33) 

< 2 , fc2 (x) ~ 4(zi)<(z2)4(z 3 ) (34) 
<3, fcl (x) := <0ri)4(* 2 )4(x 3 ) (35) 
<4, fe2 , fc3 ( x ) : = s k ni ( Xl )wf 2 (x 2 )wf 3 (x 3 ) (36) 

<5, fcl ,fc 2 (x) := Wn\( x l) w n 2 (x2)s* 3 (x 3 ) (37) 

KsMmW w k n \{x x )s k n2 {x 2 )w^nM (38) 

<iMMM^) : = <(zi)<(z 2 R 3 3 (z 3 ). (39) 

We call the functions u>™ a (x) generalized wavelets, where the index, to, rep- 
resents the smallest wavelet scale (largest ki,k2,k 3 ) appearing in the product. 
The index a indicates the values of k\ , fc 2 , &3 as well as which of the seven types 
of products appear in the basis function. 

Since the scaling function is the solution of a renormalization group equa- 
tion and all of the other basis functions are linear combinations of translated 
and scale transformed scaling functions, it follows that all of the basis functions 
have a fractal structure. This makes it difficult to represent the basis functions 
in terms of elementary functions that are smooth on a sufficiently small scale. 
Fortunately, while it is possible to compute the scaling function to any desired 
accuracy, this is never needed in any application. What is needed are the overlap 
integrals of products of arbitrary numbers of basis functions and their deriva- 
tives. It turns out that these integrals satisfy renormalization group equations 
and can all be computed exactly using the renormalization group equations and 
the normalization condition, without knowing the value of the integrand at any 
points. These equations can also be used to compute the integrals of products 
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of these functions with polynomials of arbitrary degree. Since the basis func- 
tions have compact support, and any continuous function can be approximated 
by a polynomial on a compact interval, it follows that integrals of products of 
these basis functions and continuous functions can be computed to any desired 
accuracy. Methods for computing these quantities are discussed in section five 
and the appendix. 



3 Free field conventions 

In this section we summarize our conventions for a free mass-fi scalar field and 
list some standard results that will be used in the subsequent sections. Since the 
free field operators exist and the wavelet basis expansions are exact, free fields 
provide a useful laboratory for understanding truncations as approximations. 

Free fields also play a role in modeling interacting theories. In the presence 
of both resolution and volume cutoffs, which are natural in the wavelet basis, 
all of the operator products in the Hamiltonain become finite sums of well- 
defined operators on the free-field Fock space, allowing one to formally work in 
the interaction representation. The renormalization group, discussed in section 
seven, provides one means to construct a limiting interacting theory with an 
infinite number of degrees of freedom. 

The Lagrangian density for a free scalar field of mass fj, is 

C{x) := - V*(ar) • V*(ar) - {i 2 <$>(x) 2 ). (40) 

Lagrange's equations are 

^-^(xj+^tx)^. (41) 
The generalized momentum is 



n(*) = ^jp = *(s) (42) 



and the Hamiltonian is 



= 1 J (Il(x)n(x) + V*(ar) • V*(ar) + ^ 2 <P(x) 2 )dx. (43) 

The field $(x,i) is normalized so $(x, t) and II(x, t) satisfy the canonical 
equal-time commutation relations 

[$(x,t),II(y,t)]=iJ(x-y). (44) 

The classical energy-momentum tensor is 
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d"$(x)ff'$(x) - ^r ] flv {d a ${x)d a ${x) + m 2 $(x) 2 ) (45) 
and the angular momentum tensor is 

M ^a = ^ T ua _ x » T i«xy ( 46 ) 

These tensors satisfy the conservation laws 

fyT"" = d a M^ a = (47) 

leading to the conserved charges, which after quantization, are the infinitesimal 
generators of the Poincare group 

P" = [ dxP»(x)dx V»(x) :=T^(x) (48) 

J"" := [ dxj^ v {x)dx J^(x) := M^°(x). (49) 
Jt=o 

The generators can be expressed as integrals over operator-valued densities 

H = P° = 

[ : (n(x)<b(x)-£(x)) :dx=\[: (n(:r)n(x)+V$(x)-V$(x)+^ 2 $(x) 2 ) : dx 
Jt=o 2 J 

(50) 

P = - f : n(a;)V$(x) : dx (51) 
Jt=o 

jij _ _ e ijkjk 

I dx : (n(x)x^d^(x) - n(x)x i d j ^{x)) : (52) 
Jt=o 

r° = - = K i = 

f : (\li 2 {x)x i + Jjc'VSfc) • V$(x) + l-x l fi 2 ^(x)A : dx (53) 
Jt=o \2 2 2 J 

where the :'s indicate that the operators are normal ordered. 
The field has the Fourier representation ( — h ++ metric) 

*(«) - 7^ / -^=(e-^at(p) + e -a( P) ) (54) 



\/2wp(p) 

- '/'/' 



n(x) = . ' / dp \/^4^( e " ip ' Xat (P) - eip ' x «(p)) ( 55 ) 



where 



w„(p) = vV+P 2 (56) 
is the single-particle energy. The equal-time commutation relations imply 

Kp),aV)]=<5(p-p'). (57) 
11 



The vacuum state of the held is the solution of 

o(p)|0)=0. (58) 

The infinitesimal generators have momentum-space representations as integrals 
over momentum densities 

H = Jn{p)dp = y"dpat(p)u ro (p)o(p), (59) 

P = J V(p)dp = J rfpa+(p)pa(p) ) (60) 

J = / J{v)dp = J <W(p)(i^ x p)a(p), (61) 

K = J K{p)dp = J dp a t (p )I{ i A ;CJm (p) }a (p) (62 ) 
where we set a; = in the last expression. 

4 Wavelet discretized fields 

The wavelet basis discussed in (|32M39[) can be used to write exact expansions of 
the fields as sums of products of basis functions with time-dependent operator 
coefficients 



$(x,t) = ^$ fc (n,i)^(x)+ *'KMVL(x) (63) 

n n,a,Z>fc 

n(x,i)=]Tn fe (n,^(x)+ ]T n , (n,a,t)< a (x) (64) 



n,a,l>fc 



where the operator coefficients are projections of the field operators on the 
orthonormal basis of scaling functions and wavelets 

$ fc (n,t) = f dx$(x,t)8*(x) (65) 

$ (n, a, t) = J dx$(x,t)w l na (x) (l>k) (66) 
U k (n,t) = J dxll(x,i)4(x) (67) 

tf(n,a,t) = j dxll(x,t)w l na (x) (l>k). (68) 

In these expression we only integrate over the spatial coordinates so the operator 
coefficients are time dependent. It is also possible to integrate over all space-time 
coordinates. 
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These expansions have been separated into operators smeared over scaling 
functions and generalized wavelets. This is a separation of scales - the scale-l/2 fc 
scaling functions give the coarse scale structure of the field, while the wavelet 
parts include the structure on all smaller scales. 

The field operators smeared with the scaling functions represent the average 
value of the field over a compact region weighted by the product of scaling 
functions. Because the scale l/2 fc+1 scaling function can be expressed in terms 
of the scale l/2 fc scaling functions and scale l/2 fe wavelets, (j2"4"]l . the scale l/2 fc 
wavelet and scaling-function smeared fields, contain the same information as 
the scale l/2 fe+1 scaling function smeared fields. 

When all of the wavelet contributions are included the expansions are exact. 
In any finite region there are still an infinite number of wavelet basis functions 
associated with arbitrarily small scales. 

For fields satisfying canonical equal-time commutation relations, the or- 
thonormality of the basis functions implies the equal-time canonical commu- 
tation relations for the discrete fields: 

[<S> k (n,t),$ k (m,t)} =0 [U k {n,t),U k (m,t)} = (69) 
[$ fe (n,i),n fc (m,i)] = i5nm (70) 

[$>,a,t),$ s (m,/M)] =0 [U r (n,a,t),U s (m,p,t)} =0 (71) 

(72) 

[$ r (n,a,t),$ s (m,t)] =0 [IT(n, a, t), IT(m, t)} = (73) 

[$ r (n,a,t),ll s (m,t)}=0 [lT(n,a,i),$ s (m,i)] = (74) 

where in all of these expressions k < r, s. 

Given these discrete field operators we can construct discrete creation and 
annihilation operators 

a fc (n, t) := -^=(^ k (m, t) + i^-Tl k (m, t)) (75) 

b r (n,a,t) :=^((V7<^(n,M)+^n fc (m,a,t)) (76) 
V2 

that satisfy 

[a k (n,t),a k Hm,t)]=S mn (77) 

[b\n,a,t),b>^{m,pt)] = 5 mn S 3l S a0 (78) 

with all other commutators vanishing. In these expressions 7 is a constant that 
can depend on the discrete indices. While these quantities and their adjoints 
satisfy the commutation relations (|77I78[) for creation and annihilation operators 
for any choice of 7, a specific choice is required if the annihilation operator is 
to annihilate the free field vacuum. This requires that 7 is chosen so 

(0\a k \n,t)a k (n,t)\0) = (79) 
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<0|6 rt (n,a,t)6 r (n,a,i)|0) = (80) 

(0|a fe (n,i)a fet (n,£)|0) = 1 (81) 

(0\b r {n,a,t)b r \n,a,t)\0) = 1. (82) 

With this choice of 7 the annihilation operators annihilate the mass n free-field 
vacuum. 

The solution of these equations for the scaling function fields is 



1 ± y/l - 4(0|$ fc (m, i)$ fc (m, i) |0> (0|II fc (m, f)II fc (m, f)|0) 
7( ' ~ 2(0|$ fe (m,i)$ fc (m,i)|0) 

= 1 ± y/i - 4(o|$fc(o, pgo; t)|o) (o|nfc(o, t)nfe(o^|oy 

2(0|$ fe (0,i)$ fe (0,i)|0) 1 J 

and for the wavelet fields 

_ 1± V 1 -4(0|$ r (m,Q,i)$ I '(m,a,t)|0)(0|n'-(m,Q,i)n'-(m,a,0|0) 
7( ,a) _ 2(0|$''(m,a,0$ r (m,a^)|0) 

= 1 ± y/1 - 4(01^(0, a, t)^(0, a, ^)|0) (0|n^(O ; a, t)IT(0, a, gg 

2(0|$ r (0,a,t)$f(0,cM)|CI) ' 1 ' 

These coefficients depend on the scale l/2 fe , a and the mass //; however they 
are independent of m and t by the space-time translational invariance of the 
vacuum. With this choice of 7 the scaling function fields can be expressed in 
terms of the creation and annihilation operators as 

$ fc (m,t) = -^=(a fc t( m ,t) + a fc (m,;)) (85) 
V27 



I*(m,t) = yi(o*t(m,t)-«* 
with analogous expressions for the wavelet fields, 



n fc (m,i) = i W^(a fe f(m,i) - o*(m,t)) (86) 



$ fc (m,a,i) = — L(& fct (m,a,^) + 6 fe (m,a,t)) (87) 



27 



n fe (m, t) = *y|(& fe t (m, a, t) - 6 fe (m, a, t)). 



(88) 



The coefficients 7 depend on the mass term in the field operators. This is 
because the following integrals appear in (j83H84[) : 

(0|$ fe (m,i)$ fc (m,t)|0) = -Lj / 3 n(x)4(y) eip . (x _ y)dxdydp 

(27r) d J 2w M (p) 

(0|n fe (n,i)II*(n,t)|0) = 7A3 / 5 "( x ) s "(y)^(P) e 'P-(x-y) dxdydp (90 ) 
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(0\<f r (m,a,tW(m,a,t)\0) = [ <^)K,n(y) eip . (x _ y) (gi) 

<0|n>, a ,i)IT ( n )M )|0) = / <'° (x) <" (x) ^ (p) e ^-(- y )rfxrf y rfp. 

(92) 

These are integrals of the basis functions over the two-point m&ss-fj, Wightman 
function of this field at fixed time. 

Using (|85H92I) in (|63I64[) gives exact expressions for the field operators in 
terms of the discrete creation and annihilation operators: 

^x) = ^^J=s k m (x)(a k ( m ,t) + a k H m ,t)) 

+ ]T ^-L ^wty'K «. *) + «. *)) (93) 

n(x) =<^^5i(,)(a k t(m,t) - a fe (m,t)) 

+ 4 E V^^<« (x) ^ (m ' «.*))• (94) 

In this notation the a and operators represent scale l/2 fc degrees of freedom, 
while the b and operators are associated with smaller-scale degrees of freedom. 

The Hilbert space for this free field is generated by taking limits of finite lin- 
ear combinations of products of the discrete creation operators, a fc (m, f)t, b l (m, a, 
applied to the vacuum at a fixed time. 

It is useful to decompose operators into parts that only involve the a and 

operators, parts that only involve the b and b^ operators, and mixed terms 
involving product of at least one operator from each of the above groups. 

The terms with only the a and a* operators represent the 2 _fc -scale physics, 
the terms with only the b and represent the part of the operator associated 
with scales finer that 2~ fc that do not couple with the scale 2~ k operators, and 
the mixed terms are responsible for the coupling of the 2~ fc -scale degrees of 
freedom to the smaller scale degrees of freedom. 

The Hamiltonian (|50|) has the decomposition 



H = H(a) + H(b) + H(ab) (95) 

where 

§(E ; n fc (n,0)n fc (n,0) $ fe (m, 0)^ n $ fe (n, 0) : 

n mn 

+ /i 2 ■ $fc (n, 0)$ fc (n, 0) :) (96) 

n 

HQ>) := i(£ : rf(n, a ,0)ir(n,a,0) : + 

na.i 
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53 :$'(m,a,0)D« ainij8 ^(n si 8,0): 

ma,in/3j' 

+ /i 2 ^ :&(n,a,0)&(n,a,0) :) (97) 

/,n,a 

tf(a&) := ^ E : <f'(m, a, 0)^^(11,0) : (98) 

rnrv . /n 

and the coefficients D^ n , Q n g and n are given by 

C = /dxV4(x).V4(x) (99) 

D m,a,n li9 = / ^V< a (x) • V<^(x) (100) 

D*, a>n = 2 J dxV< iQ (x) ■ V^(x). (101) 

The derivatives of the basis functions in (|99H101j) exist for Daubechies K > 3 
scaling functions. The computation of these integrals, using the renormalization 
group equations and normalization condition, is discussed in the next section 
and the appendix. These integrals are almost local in the sense that they vanish 
when the support of the basis functions do not overlap. From the above equa- 
tions we see that the part of the Hamiltonian, H(ab), that couples scale l/2 k 
to the smaller scales proceeds through the terms in the free-field Hamiltonian 
with the spatial derivatives. 

There are similar decompositions for the linear momentum (fSTj). angular 
momentum (|5"2"|) , and Lorentz boost generators (|53p . For the linear momentum 
the decomposition is: 

P = P(o) +P(6) +P(ab) (102) 

where 

P(o) = - J] : n fc (m, 0)P^ n $ fe (n, 0) : (103) 

mn 

P(6) = - J] :tf(m,a,0)P^ Qin ^(n,/3,0) : (104) 

ma.ln.l3 ,j 

P(a6) = -53 itf^a^JP^/^O):, (105) 

ma,/n 

and 

PL,n = -/HWV4x) (106) 

P»,a,-,fl = - / ^L, a (x)Vx< p (x)+ (107) 
PL.a.n = - J ^(xJVx^x) + »;(x)V I tt/ n)0 (x)). (108) 
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For the angular momentum 

J = J (a) + J(b) + J(ab) 

where 

J(a)=£:n*(m,0)J mn <I> fe (n,0): 

mn 

J(6)= Yl :n'(m,a,0)J^ a;n;/3 ^(n,/3,0): 

ma, in/3 ,j 

3(ab) = £ ■■ tf(m, a ,0)J^ Q!n $ fc (n,0) : +II fc (n,0)j£ mia $ , (m> 

ma Jn 

and 

J m n = - J dx( Sm (x)(x x V x ) S £(x) 

J m,a,„,/3 =~ J d™L, a W(X X V X K^(X), 

JL,„ = - J dx< !0 W(x x V x )**(x), 

JL,„ = - y*dxa*(x)(x x V x )w; m> a(x). 
The decomposition for the rotationless boost generators is 

K = K(a)+K(fe) + K(a6). 

where 

K ( fl )= / ^E( F nmOn(n,0)n(m,0): 

•* m,n 

+/i 2 : $(n,0)$(m,0) : +G nm : $(n,0)$(m,0) :) 
K ^=jl £ ( F m, Q ^0n'(m, a )ff(n,/3,O): 

+ M 2 : ^(m,a,0)^(n,/3,0) : +G^ an/3 : $'(m> a, 0)<F(n, /?, 0) :) 
K(a6)=/ £ (F ma , n (:rf(m,a)n fe (n,0): 

Z , m , a , n 

+/i 2 : *'(m,a,0)$ fc (n ) 0) : +G£ Q>m : *'(n, a, 0)$ fe (m, 0) :) 
and the overlap integrals are 

FL - / dxs*(x)x£(x) 
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G^ m ^|dxV4(x)-V4(x)x (122) 

F l L mf 3 = / dx<Jx)x<,(x) (123) 

G l L mf 3 = J dxVw l aa (x) • V«4g(x)x (124) 

m = J dx< a (x)x&(x) (125) 

GLL = | dxV< a (x) • V&(x)x. (126) 



All of the numerical coefficients can all be computed exactly using linear algebra 
and the scaling equation. They vanish when the supports of the wavelets or 
scaling functions have no overlap. 

5 Integrals 

The expressions for the Poincare generators derived in the previous section 
are linear combinations of normal products of discrete wavelet fields multiplied 
by numerical coefficients. The numerical coefficients are integrals of products 
of scaling functions, wavelets, first derivatives of these quantities and powers 
of x. For a 3 + 1 dimensional field theory each coefficient is a product of 
three coefficients involving one-dimensional integrals. In this section we discuss 
methods for computing these one-dimensional integrals. 

While all of the integrals involve functions with fractal values, making stan- 
dard numerical methods impractical, the unique properties of the wavelet basis 
makes it possible to compute all of the integrals appearing in the previous sec- 
tion exactly. In this section we discuss how to compute the coefficients in the 
previous section. We follow the general discussion with an illustrative example. 

The integrals in the last section are products of one-dimensional integrals of 
the form 

T„= J h(x)f 2 (x)---f n (x)dx. (127) 

where the functions fi(x) are scaling functions, wavelets, first derivatives of 
these functions or powers of x 

The important observation is that all of the functions that appear in the 
integrand of the coefficient r„ have simple transformation properties under scale 
transformation and translations. This, along with the scaling equation and 
normalization condition, can be used to reduce the computation of all of these 
quantities to finite linear algebra. 

The computation uses the following relations, which follow from the defini- 
tions: 




(128) 
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DT = T D (129) 

± D = 2d4- (130) 
dx dx 

Dx = 2-SxD (131) 

Tx = (x- 1)T. (132) 

In addition, we use the scaling equation, the definition of the wavelet and the 
derivatives of these equations in the form (|25H27I) : 

s k m {x)=Y J H mn s k n +1 (x) (133) 

n 

■ n (x)=J2 G ^n +1 (x) (134) 

n 

2^iWl fc+1) ' (135) 

n 

:=2^G mnS ( fc+1) '- (136) 



k 

w 



These equations can be used to express s^x), w^x), sjj(a:), and w k [(x), with 
scale 2~ fc , as linear combinations of scaling functions and their derivatives with 
scale 2~( k+1 \ Repeated application of these equations can be used to increase 
k in each of the functions by any desired amount. 

In addition, the scale factor k of all functions in the integral can be increased 
or decreased by the same amount using 

j D k f 1 (x)D k f 2 (x) ■ ■ ■ D k f m (x)dx = 
2 Mf -i) f f 1 ( x )f 2 (x) ■ ■ ■ f m (x)dx (137) 



where in this equation k can be positive or negative. We recall from the defi- 
nition Q and ([S]) that D increases k by 1 when applied to scaling functions or 
wavelets. For derivatives we also need to use (|130[) . 

Using (|137j) to make all of the fc-values negative and then using (|133H136j) to 
increase each one to zero leads to an expression for the integral (|127j) as a linear 
combination of integrals involving products of scale 1/2° scaling functions, their 
derivatives and powers of x. 

It is useful to first consider integrals where there are no powers of x; we will 
see later that the corresponding integrals with powers of x can be expressed in 
terms of the integrals with no powers of x. To compute a general coefficient, r n 
we use the following steps. 

1. Step 1: Use (|1 3T[) to relate the integral to another integral where the finest 
scale appearing in the integrand is 1/2°. Then the scale of each function 
in the integrand is 1/2° or coarser (negative k). 
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2. Step 2: Use repeated application (I133fll36|) to replace all of the coarse 
scale functions by linear combinations of scale 1/2° functions. The result 
of these two steps is that original integral can be expressed as a finite sum 
of coefficients where all of the functions in the integrand are scale 1/2° 
scaling functions or their derivatives. Note that the wavelet contributions 
can always be expressed in terms of scaling functions using (|134l or H36p . 

3. Step 3: Use integer translational invariance to make the support of the left 
most function start at 0. Then, because each function has compact support 
on an interval of width [2K-1], non-zero coefficients have translational 
indices that vary from [— 2K + 2,2K — 2]. That means the number of 
non-zero coefficients that cannot be generated by translations is less than 
(4JC- 3)"- 1 (finite). 

4. Step 4: Apply the scaling equation to the integrand of the resulting inte- 
gral and use translational invariance to shift the left most index to zero. 
This results in a set of (AK — 3) n_1 homogeneous linear equations for the 
coefficient functions with left index zero. 

5. Step 5: Use the partition of unity property (there is also one for the deriva- 
tive of the scaling function) to get additional inhomogeneous equations for 
these coefficients. 

6. Step 6: Combine the homogeneous and inhomogeneous equations to get a 
linear system and solve for the coefficients. 

7. Step 7: Reverse the steps, using the finite number of solution of the linear 
equations to construct the general coefficients. 

All of the steps are straightforward; for step six some of the equations are 
redundant and one must select a set of independent equations that include at 
least one inhomogeneous equation. 

Wc illustrate construction with the example of computing the integral of the 
product of a scaling function, wavelet, and derivative of a scaling function with 
different scales 

f ds k+3 
I k rtf' k+j ■= J s k m (x)w k n +l (x)^ r (x)dx. (138) 

To be specific we assume that j > I and k > 0. For the first step we use (|137[) 
and lfT3"0"f to reduce the finest scale, l/2 k+j to 1/2° 

jk,k+l,k+j _ 2 -3(fc+i)/2 r- j,-j+l,0 fion\ 
mnp mnp ' \ '-'^J 

This makes the last scale index and the other two negative. Since all of the 
other scale factors have a negative exponent, we can use (1134II135| to increase 
each of them to zero: 

T-j,-j+l,0 _ \ " ttJ ( G m-!-h rO.0,0 / 14Q n 
1 mnp — / j 11 mm'\ { ~ r11 M»' I m , n'p K 1 ^) 

m'n' 
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where in (|140l) H ] mm , is the product of j factors of the matrix H mn . While 
the matrix H mn is formally infinite, only a finite number of terms in this sum 
contribute to a particular values of m, so the sums in (|140[) are finite. Similar 
remarks apply to G-fP _1 . The integral on the right side of (|140[) is 

Cn' P = J S m ,{x) Sn ,(x)^{x-p)dx. (141) 

Note that the wavelet in (|138H139j) is replaced by a scaling function because 
equation (|134|) expresses the scale 1/2' wavelet as a linear combination of scale 
l/2 i+1 scaling functions. 

Next we use translational invariance to transform the first index to zero 

jo.o.o =J o,o,o (U2) 

m'n'p 0,n' — m' ,p— m' ' V J 

Thus we can express any of the integrals, Imnp l ' k+:> ' > m terms of the (2K — 2) 2 
non-trivial coefficients / mn • 

Because of the compact support of the scaling function and its derivative, 
^Omn i s non- vanishing only for values of m, n satisfying 

- 2K + 1< m, n < 2K + 1; |m - n| < 2K - 1. (143) 

These integrals have the form 

lamn = J s(x)s(x - m )^( x ~ n)dx. (144) 

Next we use (|133|134|135|) and (|137j) and translational invariance to get the 
following scaling equation 

r 0,0,0 _ 1 J- 1,-1,-1 _ 2 \^ TT TT TT fOOO _ 

-'Omn ~ /7; J 0mii ~ /77 / , ' 12 mm ' n nn' * l> m ' ' n' ~ 



I m n 

V2 ^2 H u , H mm > H nn , Io° r °, _ t , n ,_ v (145) 



I'm'n' 



which is a system of (2K — 2) x (2X — 2) homogeneous equations for the values 
of 1^ with -2K + 1 < m, n < 2K + 1. 

In order to solve this system we also need an inhomogeneous equation, which 
must be related to the normalization condition on the scaling function. The 
Daubechies K — 3 scaling function is defined so 1, x and x 2 can be pointwise 
expanded as linear combinations of scaling functions 

1 = y^^a n s n (x) and a; = ^ 6 n s„(x). (146) 

n n 

The normalization condition ([3]) and the orthogonality condition (|16p give a n = 
1 or 

l = J2sn(x), (147) 
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which is a restatement of the partition of unity property. To calculate b n use 
the unitarity of D and the scaling equation to get 

xs n (x)dx = n + / xs(x)dx = n + / D~ 1 xD~ 1 s(x)dx = 



n-\ — hi f xs(x — l)dx — 

2%/2 ^ J 

!i+ ^~f^Yl lhl J s ( x ) dx + 7^J2 hl J xs ( x ) da 



u —j= Tjhi + -J xs(x)dx (148) 



2^2 , 

where we used CL2]), J^i h i = >/2, in ([US]) . Setting n = gives 

/ xs(x)dx = —= Zft./. 
7 V 2 j 



(149) 



It follows that 



and 



6„ = n + 



^=J2 lhi ( 15 °) 



x = y^ j ns n {x) + y^Z/i;. (151) 

If we differentiate (II 5 1 p we get a partition of unity for the derivative of the 
scaling function, 

l = Y j ns' n {x). (152) 
Equation (|152j) implies that the coefficient functions satisfy 

n 

This gives a non-trivial inhomogeneous equation when m = 0. If we sum over 
the first index we get the equivalent inhomogeneous equation 

^n/omn = l- (154) 

The next step is to identify a set of independent linear equations, including 
(|154p , for the non-zero InLL. Solving these equations gives the non- vanishing 

r0,0,0 
Omn ' 

Reversing these steps we get the following expression for the general coeffi- 
cient in terms of the coefficients ij^j : 

T k,k+l,k+j _ 2 _3(fc+j)/2 >p H j (nfrj-i-i) ,7-0,0,0 
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This shows that the scaling properties lead to explicit formulas for an infinite 
number of coefficients in terms of the finite set of coefficients /q m°n with —2K + 
1 < m, n, m — n < 2K — 1 . 

Given the above method to compute all of the integrals without powers of x, 
is it easy to generalize the construction to treat integrands involving powers of 
x. This computation uses the translations that appear in the scaling equations 
to reduce the desired integral to one of the types discussed above. 

The basic trick is to combine (jl33H137[) . For example, all of the moments, 



< x m > s = / s(x)x m dx < x m > w = J w(x)x m dx, (156) 
can be constructed recursively from the normalization condition 

< a; > s = / dxs(x) = 1 (157) 



using the scaling equation 

< x rn > s = I D- l x m D- l s{x)dx 



1 1 

71 2™ 



J2 h l f {x + l) m s{x)dx 

= 4=^- y h Y — m! ,,, r- fc / x k s{x)dx. (158) 

Using hi — y/2, and moving the k — m term to the left side of equation (|158|) 
gives the recursion relation for the m-th moments in terms of k < m moments 



m-l , fiK-X 



J xms ^ dx = 2^171 S ( £ hll j I xks[x)dx - (159) 

Moments of the wavelets are obtained by replacing hi in (|158[) by gf. 



< x rn > w = J w(x)x m dx = 

= 9lY Tu m! ,x, ?T "" fc / X k s(x)dx. (160) 

s f2 2 m ^ y k\(m- k)\ J 

This means the first moment can be expressed in terms of the th moment, which 
is 1 by the scale fixing condition. This method can be generalized to calculate 
moments multiplied by arbitrary products of scaling functions, wavelets and 
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their first derivatives. In the general case one gets a linear system for the k in 
terms of the 0, 1, • • • , k — 1 st moments. 

Using these methods all of the coefficients appearing in the expressions for 
the Poincare generators can be expressed in terms of the scaling coefficients, hi, 
and the solution of a finite system of linear equations. The factors of £P and 
GiP - ' -1 in (| 1 55|) are dominated by the terms involving j or j — I — 1 products 
of hi = g + 4 « 1.1 or H mn -> l.l<5i_ n _ 2TO and G mn -> l-1^4_ n _ 2m . 

Calculations of (|99H10ip are done explicitly in the appendix. 

6 Poincare Symmetry 

There are three properties of the wavelet basis that are useful for investigating 
Poincare symmetry breaking in discrete truncations of field theory. First, the 
scaling functions on a fixed scale are a compactly supported, locally finite par- 
tition of unity. This makes it possible to test the Poincare symmetry in finite 
volumes. Second, the scaling functions and wavelets form a basis, which means 
that it is possible to precisely identify the corrections that are needed to restore 
Poincare invariance in truncated theories. Third, because it is possible to com- 
pute the coefficients of the operators that restore the commutation relations to 
a set of truncated generators, it is possible to identify and classify the size of all 
of the correction terms. 

In this section we discuss these three properties in more detail. First we 
consider the partition of unity property. Up to an overall multiplicative constant 
the scaling functions on a given scale, l/2 fe , form a partition of unity that has 
the form 

1 = ^2-3*/^). (161) 
m 

While this sum is formally infinite, at any given point there are at most (2 K— l) 3 
non-vanishing terms in this sum. 

The partition of unity can be utilized to decompose each of the Poincare 
generators into sums of localized operators. To do this note that Noether's 
theorem gives formal expressions for the Poincare generators as integrals of local 
densities (|48H49[) over the t — surface. The local densities involve products of 
field operators and their derivatives at the same spacetime point. 

The local density, O l (x), of the i-th Poincare generator satisfies 

[O i (x,0),OJ'(y,0)]=^(x-y)/« fc O fe (y ) 0) (162) 

where the f^ k are the structure constants of the Poincare group. Generators 
are integrals of these densities over the t = hyperplanc 

Oi := J dxO l (x,0). (163) 

These integrals are the conserved charges in Noether's theorem and formally 
become time independent after integrating over a fixed-time hyperplanc. The 
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Poincare commutation relations are obtained by integrating both sides of equa- 
tion (|162p over a fixed-time hyperplane. This gives the commutation relations 



[O i ,O j ]=if ijk O k O l = J dxO l {x,0). (164) 

While the generators are not local operators, they involve integrals of local 
densities over a 3-dimensional hyperplane, and the partition of unity for scale- 
l/2 k scaling functions can be used to express the commutation relations in an 
almost local form. To do this we insert the partition of unity, (|161[) . in (|163l) to 
get 

0* = J2 I °m (165) 

m J 

where the smeared density is 

O m ■■= J 2- 3fc / 2 4 1 (x)O i (x,0)dx (166) 

and the integrals are over the t = hyperplane. 

In this notation the exact commutation relations (|164[) can be expressed in 
terms of the sums: 

£[0* 0£']= O*'. (167) 

mn I 1 1 

While these sums have an infinite number of terms, there are no convergence 
problems because when (|167[) operates on the dense set of states |^) constructed 
out of finite linear combinations of products of creation operators smeared 
against wavelets and scaling functions, whose support is necessarily contained 
in a bounded region, then only a finite number of terms in the sum (|167|) are 
non-zero, giving 

mn I n 

In (|168[) the sums are finite, with m and n running only over the indices corre- 
sponding to scaling functions in the partition of unity that have common sup- 
port with the smearing functions s n (x) and iuj(x) used in the discrete creation 
operators. 

Equation (I168p is an exact consequence of Poincare invariance, and if it is 
satisfied on the dense set of vectors generated by finite numbers of discrete 
creation operators, one recovers exact Poincare invariance. 

Thus Poincare invariance can be tested by replacing the infinite sums in (|167l) 
by a suitable finite sum of the localized operators, O n '. Where approximations 
enter is that the local densities 0*(x, 0) that appear in the definition of O^ 1 are 
constructed out of normal products of field operators. The fields that appear in 
these expressions can all be expanded in the form (|63H64[) . The partition of unity 
projects the product of different scale discrete fields on the scale l/2 fc subspace. 
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As long as there is no truncation everything is exact. The problem is that for 
normal products of more than one field, there are scale l/2 fe contributions to the 
operator products that involve fields on scales smaller than l/2 fe . If we truncate 
the theory by discarding all of contributions to the field operators associated 
with scales smaller than l/2 fe , then the discarded scale l/2 fc contributions that 
come from projection of products of smaller scale field operators are precisely 
the corrections that are needed to restore Poincare invariance. 

To understand the structure of these discarded terms consider the exact 
expression for the linear momentum of a free field using the scale l/2 fe scaling 
functions and associated wavelets. The scalar density 0(x, 0) for the momentum 
operator is (fSTj) which is the normal product the n(a;) field and the gradient of 
the $(x) field: 

P(x) = - : II(x, 0)V$(x, 0) : . (169) 
Using the partition of unity (|161l) equation ([5"Tj) can be expressed exactly as 

P = E P m (170) 

111 

where 

Pm = " / ■ H(x, 0)V<£(x, 0) : 2- 3k ^ S k m (x)d X . (171) 

Alternatively, we can express using the exact expansion of each of the 
fields in the scale l/2 fe wavelet basis and smearing the result with the scale l/2 k 
partition of unity: 

Pi=E:n fc (j,0)P^$ fc (n,0): + 
E :n'(j, a ,0)^ fe Q fe nm $ fe (n,0): + 

i,Q,j,n 

E :n fc (n,0)^ fe Qjm <i>'(n,a,0): + 

!,a,j,n 

E :n ; (n )a) 0)P/tn, /3 , m ^ s (n,A0): (172) 

/,a,j,n 

where the expansion coefficients are 

P%& = -/ 4(x)^(x)V Sj fe (x)dx (173) 

Pt~ = - J t4 a (x)«*(x)V£(x)dx (174) 
P?!n k a , m = ~J S L(x) S j s (x)V<, Q (x)dx (175) 
/tL.v.,,, = - J ^ iQ (x) S i(x)V<^(x)dx. (176) 
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The coefficients above are numerical coefficients. 

To understand the interpretation of these results we note that f| 168[) is an 
exact consequence of the commutation relations expressed in terms of the scale 
l/2 fc partition of unity. The expansion (|172l) is exact, and the resulting operators 
satisfy (I168[) . If the fields are replaced by resolution the l/2 fe approximations 
given by first term on the right of and (IM|) , then (II 72[) would be replaced by 
an "approximation" consisting of only the first line of (|172l) . This approximation 
would violate the commutation relations (I168p . 

The "corrections" defined by the second through fourth line of (|172p provide 
the missing physics from smaller scales that is needed to restore the commu- 
tation relations. Note that (|168l) only involves operators that act in a finite 
volume, so any problems due to a volume truncation are not relevant. 

In the free field case the generators (|59H62p are sums of normal products 
of one creation and one annihilation operator. The commutator of two oper- 
ators with the structure a k (m) ji a k (n) leads to another operator of the same 
form. In addition commutators of operators of the structure a fe (m)'l'6'(n, a) 
with b l (n, a)^a k (j) also lead to structures of the form a k (mya k (j). 

The breaking of Poincare invariance is related to the question of how large 
are the corrections P^ nm , P^aa m an d ^janflm that restore the commutation 
relations relative to the coefficients PjJJj . that define the approximate genera- 
tors in a model with a resolution l/2 k cutoff. 

These coefficients can all be computed exactly using the methods of the 
previous section in order to identify the largest correction terms. 

7 Renormalization Group 

The wavelet basis decomposes the Hilbert space L 2 (]R 3 ) into an orthogonal direct 
sum of infinite dimensional subspaces of successively finer resolution. 

The Hamiltonian can be truncated at any scale. The resulting operator 
is something like averaging over lattice blocks, except the block functions are 
replaced by product of scaling functions of a given scale. In our free field example 
the scale l/2 fe truncated Hamiltonian has the form H(a) in (|96l) . 

It is useful to consider the relation of the free Hamiltonians H (a) truncated 
on two different scales. To do this we use the scaling equation in the form (|22[) 
to obtain the following identities: 



6; 




s' 



k+i 





(177) 



and 
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h-imh^n f s k+v (x)s k+v (x)dx = ht^mhj^nDfr 1 . (178) 
jl jl 

Next consider the individual terms in the resolution l/2 fc+1 truncated free 
field Hamiltonain H k+1 (a). This Hamiltonain is the sum of the following three 
terms 

i5^:n*+ 1 (n J 0)n fc + 1 (n > 0):, (179) 

n 

■■ $ fc+1 (m,0)^+ n 1 $ fe+1 (n,0) :, (180) 

mn 

and 

fi 2 J2 ■ ^ +1 (n,0)$ fc+1 (n,0) :). (181) 

n 

The discrete fields in these terms have the form 

$ fe+1 (n,0) = J ^0) S k n + 1 (x) S k + 1 (y)s k n t 1 (x)d X (182) 

and 

II fc+1 (n,0) = J U(^0)s k + 1 (x) S k n t 1 (y) S ^ 1 (x)d x (183) 
Using (IM1) these can be expressed in the form 

n fc+1 (n,o) = 

n(x,0)«* l (a;)«^(y)«* | (a:)dx = 
/n(x,o)(^ hnx—2]! s j 1 (x 

jl jl 
32 ji js js 

(184) 

$ fe+1 (n,0) = 
J $(x,0)4(x)4(y)4( S )dx = 

J *(X,0)(5^ h ni ^ 2jl S k 1 {x)+Y9n 1 -2 J1 W k - 1 {x))x 

jl jl 

(^2 h n2-2j 2 S k h (y)+Y 9n2-2j 2 ^ 2 _1 (2/))(5Z h n 3 -2j 3 S k 3 (z)+Y 9n 3 -2j 3 W k j~ l (z))dx. 
h h ji js 

(185) 

These expressions have the form 

$ fe+1 (n,0) = J2H hn *-^ k &V + ■■■ (186) 



28 



n fe+1 (n, 0) = H ^- 2ji n fe (j, 0) + • • • (187) 
j * 

where the • • • terms represent terms where at least one scale k wavelet w k (x) 
appears in the integral. Using these expressions in (jlT9[) gives 

i(^:n fc+1 (n,0)n fe+1 (n,0) := 

n 

: n^-^ nfe (j' o )n^- 2 ^ nfc ( r '°) +••• = 

n / m 

i(^:n fc (j,0)n fe (j,0) + ..-= (188) 

n 

where we have used (|177|) three times. We see that this term in the resolution 
l/2 k+1 Hamiltonian can be expressed as the corresponding term in the resolution 
l/2 fc Hamiltonian plus wavelet related corrections that restore the full l/2 fe+1 
scale physics. 

Similarly we find the mass term has the same form 

/i 2 £ : $ fc+1 (n,0)$ fe+1 (n,0) : ) = M 2 ^ : $ fc (n, 0)$ fc (n, 0) :+•••. (189) 

n n 

Finally we consider the term (|180p 

5^:$ fc + 1 (m,0)Z?S B 1 ^ 1 (ii J 0) := 

mn 

£:(£n^-%**tf»°) +•••)* 

mn j i 

£ ^(E II h ni -2i^ k (l, 0) + ■ • • ) (190) 

1 i 

Using (fTTH) and (IT78| this becomes 

£:$ fc (m,0)P* m **(n,0):+-.. . (191) 

mn 

Taken together these calculations show that the resolution l/2 fc+1 free Hamil- 
tonain is equal to the resolution l/2 k free Hamiltonian plus a correction that 
fills in the missing degrees of freedom that appear in the scale l/2 k+1 but not 
on scale l/2 k . The correction terms all involve fields smeared with wavelets. 
An important observation is that for free fields the derivative terms provide the 
coupling between the two scales. A similar analysis can be used to show that 
all of the free-field Poincare generators have the property that the resolution 
l/2 fc+1 operators are equal to the resolution l/2 fe generators plus wavelet field 
corrections. 
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While we derived these results by considering the example of a free field, 
the result also holds for interactions. For example a scale l/2 k+1 truncated 
J : $™(x, 0) : dx interaction has the form 

:$ fe +Vi,0)---$ fe+ Vn,0):r^+ 1 .. mii (192) 

tni ■ ■ -win 

where 

rL + 1 1 .. m „ = J 4 + / (x) ■ • • (x)dx. (193) 
It is straightforward to show, using (|2"2"]l and (|2~i]) . that 

r mi- - m„ = H h hi-2 mii ■ ■ ■ ^„«-2m„ i rf ] + . 1 lti (194) 

which when used with (jl86[) and (|187l) gives 

£ : $ fe+1 ( mi , 0) ■ • • $ fc+1 (m n , 0)lt+!.. mn = 

mi ■ ■ -TOn 

^ :* fc (m lj 0)-..$ fc (m n ,0)rS 11 ... niB +... (195) 

mi ■ ••m n 

where the • • • terms represent the contributions where at least one of the ex- 
pansions functions is a scale l/2 fc wavelet, w k n (x). 

It follows that a resolution l/2 fc+1 Hamiltonian has the form 

H k+1 {a) = H k {a) + H k (b) + H k {ab) (196) 

where H k (a) is the part of H k+1 that has only scale l/2 fc scaling function fields, 
H k (b) has only scale l/2 fe wavelet fields and H k (ab) contains the terms with 
at least one wavelet and one scaling function field. In the interacting case the 
interaction also contributes to H k (ab) 

The creation operators associated with scale l/2 fe scaling function fields and 
the creation operators associated with at least one scale l/2 fc wavelet field each 
generate mutually orthogonal subspaces on the Fock space generated by the 
scale l/2 fc+1 scaling function creation operators. 

On can then seek a unitary transformation that Block diagonalizes H k+1 (a) 
on the subspace generated by the scale l/2 fc scaling function creation opera- 
tors. The resulting effective Hamiltonain will be a Hamiltonian in the scale 
l/2 k degrees of freedom that includes the effects of the eliminated scale l/2 fe+1 
degrees of freedom. Both Hamiltonians are different functions of the same cou- 
pling constants and mass parameters. Renormalization of the parameters in 
the Hamiltonian of the model is necessary to keep some resolution l/2 fc observ- 
ables fixed. This gives a new coarse-scale Hamiltonain that includes additional 
degrees of freedom at a finer scale. 

Because the resolution l/2 fc Hamiltonain has the same form for any k, the 
starting scale is arbitrary so we can repeat this process, successively eliminating 
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degrees of freedom associated with smaller and smaller scales, renormalizing the 
parameters in the Hamiltonain at each step. At some point one can stop and 
the result will be an effective theory that describes the scale l/2 fc degrees of 
freedom, including the effects of the eliminated scale l/2 k+l degrees of freedom, 
or one can proceed to try to find a fixed point of this renormalization group 
equation. 

The are a number of possible approaches that can be used to eliminate the 
scale l/2 fe+/ degrees that appear in the scale l/2 fe+1 Hamiltonain but not in the 
scale l/2 fc Hamiltonian. 

An efficient method is to use the similarity renormalization group method [TS]. 
This involves solving the differential equation 

^& = [H(X)[H(X,H k (ab)}} (197) 

with initial condition 

H(0) = H k+l (a) = H k (a) + H k {b) + H k (ab). (198) 

The resulting H(X) will evolve to a Hamiltonian that is approximately block 
diagonal on the space generated by scale 1 /2 k scaling function creation operators 
and its orthogonal complement. In principle this can also be applied to the full 
Hamiltonian. It is easy to see the iterative solution of this equation generates 
more complicated interactions in the Hamiltonian with each iteration. The 
iteration involves commutators of discrete creation and annihilation operators 
multiplied by algebraically computable coefficients. 

A second approach is to try to do this perturbatively. In this case the 
Hamiltonian is expressed in the form 



H(a) H 1 
(H l y H c {a) 



(199) 



where in this expression H (a) is the part of Hamiltonian that maps the Fock 
space generated by the scale l/2 fc scaling function creation operators, H c is 
the projection of the resolution l/2 fe+1 Hamiltonian on the orthogonal comple- 
ment of this space, and H 1 and {H 1 )^ are the parts of the resolution l/2 fc+1 
Hamiltonian that couple these two spaces. 

A unitary transformation of the Okobu [16] [17] form 



TI _I (/. + AU)"Va -At(/ i + AAt)-Va . 



with 

a -{am o) ^ 

will block diagonailze the Hamiltonian on the Fock space generated by the scale 
l/2 fc creation operators provided A satisfies 

AH (a) - H c (a)A + (H 1 )^ - AH 1 A = 0. (202) 
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Equation (|202[) is a non-linear equation for A that can be solved perturbatively 
in H 1 . In this case we seek a solution for A = 1 of 



AH(a) -H c (a)A + \(H 



J ) t - XAH'A = 0. 



(203) 



with 



oo 



A = A " A . 



n 



(204) 



n=i 



An important feature of the wavelet method is that the coupling of the scales 
proceeds through the terms H(ab). These terms involve discrete operators mul- 
tiplied by coefficients that involve integrals over product of scaling functions, 
wavelets and their derivatives. Since these coefficients can all be computed ex- 
actly (using finite linear algebra) it is possible to identify the largest or most 
important terms and discard the smaller ones to get more efficient approximate 
solutions. 

8 Gauge Invariance 

To motivate the implementation of local gauge invariance in the wavelet rep- 
resentation of field theory we consider the example of an SU(3) gauge field. 
The treatment of full gauge invariance with respect to fields smeared over four 
space-time variables or time-independent gauge transformation restricted to a 
fixed-time hyperplane for fields smeared over a basis for the hyperplane is sim- 
ilar. Since most of our development has been for fields smeared over a basis for 
the fixed-time hyperplane, we limit our discussion of gauge transformations to 
time-independent gauge transformations. 

In the wavelet basis the basis functions are not local, however if all of the 
basis functions are retained linear combinations can be used to describe ob- 
servables associated with arbitrarily small regions. Locally independent gauge 
transformations can be built out of independent gauge transformation associated 
with each of the individual discrete wavelet fields. To be specific we consider an 
SU (3) color gauge group. 

We consider transformations of the form 



where V* c ,(n,t) and V c z c ,(n, a, t) are SU(3) valued functions of the various pa- 
rameters. They are independent for each independent discrete field operator. 
It is obvious that quadratic expressions of the form: 



<&*(n,t) *c(n,t) = £V&(n,t)**(n,t) 



(205) 



C 



$i(n, a, t) -> $^(n, a, t) ^ V l cd (n, a, t)$J, (n, a, t) 



(206) 



]T $*(n, i)$*(n, t) + *c(n, a, t)¥ c (n, a, t) 



(207) 
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are invariant with respect to the gauge transformations (|205l - 12061) . 

The construction of a representation for the covariant derivative can be de- 
duced from the continuum covariant derivative 

- VD cd = -V5 cd + A d (x, t)X d cc , (208) 

where X^ c , are the Gell-Mann matrices. Matrix elements of these operators in 
the wavelet basis have the form 

" C;cc' = -^ k mn;cC + E A mn;d(x, t)>& (209) 

d 

_ '^mnaic' = _ ^rano;cc' ^~ ^ ] ^iima;il( x i ^)^cc' (210) 

d 

~ ^man;cc' = _ ^mon;cc' ~l~ ^ ] ^mtm;(l( X i ^)^cc' (2H) 

~ ^LnftcC = -•'VLnPicC + E A - W (^ ^ (212) 



where the multi-index quantities appearing in (|209H212"1) are the matrix elements 
in the wavelet basis 



V^ n;cc , = J ^(x)V^4(x)dx (213) 

A^ n;d = W 4(x)A d (x)4(x)dx (214) 

V^ nQ;cc , = / 4(x)V5 cc ,< a (x)rfx (215) 

A^ nQ;d = E / 4WA rf (x)< a ( X )dx (216) 
i J 

V£ an;cc , = | < a (x)V<Ws*(x)dx (217) 

A^ Qn;d = E / < Q (x)A d (x)4(x)dx (218) 

Vi nftc , = J < a (x)V5 cc ^^(x)dx (219) 

A^n^d = E/ <aWA d (x)^(x)dx (220) 



This derivative will transform covariantly 



V k c ,(m,t)V k cc , =V k >V c k c ,(n,t) (221) 
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V c k c , (m, t)V k ^ cc , = 2>«' nQ;cc , Vi, (n,a,t) (222) 

V l cc , (m, a, t)V l ^ an:cc , = V l £ aIl . cc ,V k c , (n, a, t) (223) 

(n, a, tpi nj3;cc , = 0^' anftcc ,Kc' (n, A t) (224) 
provided the vector potential matrix transforms like 

-iV&(m, t)V k mn + £ Kc"(m, t)A k ma . d \ d c »c = 

dc" 

-iV k ma V k An,t)+^A' k a .^VXM (225) 

dc" 

-iV k , (m, t)V^ nQ + ^ V* , (m, t) A^ nQ;d Af„ c = 

dc" 

-iV^ na Kc'(n,a,t)+^A^ n;d Al,,K"c(n,a*) (226) 
-iV^ (m, a, i)V^ Qn + V l" ( m ' a ' *) A L;^f» c = 

dc" 

- iV£ an V£, (n, i) + £ A'^ n;d A^K"c(n, *) (227) 

( m > «. + £ K'c" (m. «. t) A Ln0; d ^c = 

dc" 

- *V£ an/3 ^< (n, M + £ A m Q „ ft dAl" K-c(n, ft t) (228) 

dc" 

With these transformation properties the quantities 

-^$ c fe (m,^ n;cc ,^(n,i) (229) 

inn 

- * ^ (m, a, t)2>£ an;cc ,$*(n, i) (230) 

man 

- i *c (m, t)2^ na;cc , # c (n, a, i) (231) 

inn 

- i ^ $i(m, a, *)©LnflicC*cK ft *) (232) 

mn 

are invariant. 

Because these are invariant, matrix element by matrix element, the invari- 
ance is preserved by truncation. A covariant field strength tensor is obtained by 
taking the commutator of the covariant derivatives projected on different axes. 

Note that the generalization to 3+1 dimensions is a direct extension of the 
three-dimensional results. The interesting feature is that in the full theory the 
gauge invariance is implemented by an infinite number of independent non-local 
gauge transformations that act independently on each degree of freedom. 
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9 Summary and Conclusion 



In this paper we discussed some of the advantages of using the basis generated 
by Daubechies scaling functions and wavelets to formulate exact discretizations 
of local field theories. We emphasized on the special properties of the basis and 
how these could be useful in field theoretic applications. The most important 
properties are that the basis functions have compact support, contain locally 
finite partitions of unity, and are related to fixed points of a renormalization 
group equation. 

The working assumption is that it is possible to make sense out of quantum 
fields smeared with a class of test functions that are not infinitely differentiable, 
but instead have a fractal character with a limited amount of smoothness. The 
justification for this is that the integral of a product of Daubechies K = 3 
scaling functions over the Kallen-Lehmann representation of a mass fj, two-point 
Wightman function exists. This justifies the use of wavelet smeared fields in free- 
field theories and more generally in models with volume and resolution cutoffs. 

Given this assumption the wavelet basis leads to an exact representation of 
the local field as an infinite linear combination of smeared fields multiplied by 
compactly supported basis functions. The smeared field operators are operators 
rather than operator valued distributions. They generate a local algebra in the 
sense that in any open set of spacetime there are operators associated with 
smearing functions that have support entirely in that open set. 

In the wavelet representation products of local field operators are replaced 
by infinite sums of well-defined operators multiplied by products of wavelet basis 
functions at different space-time points. In this representation singularities in 
the operator products at nearby points are replaced by convergence questions. 
For example, the local composite fields in the Wilson-Zimmermann formulation 
of the operator product expansion [T|] are recursively constructed by identifying 
and ordering the most singular matrix elements as the separation between points 
vanish. In the wavelet representation these matrix elements are represented 
by infinite sums of products of basis functions at different points with well- 
defined expectation values of discrete field operators. The identification of the 
singularity class of a given matrix elements depends on the asymptotic properties 
of the series as a function of the separation between points. 

The discussions in this paper focused on fields smeared only over wavelet 
basis functions in three dimensions at a fixed time. This representation is more 
convenient for dealing with the Hamiltonian formulation of the dynamics as 
well as the treatment of the Poincare symmetry. In this representation, if the 
basis is truncated to a finite number of degrees of freedom (equivalent to a 
volume and resolution cutoff), then the interactions in the Hamiltonain become 
well-defined operators and one can in principle solve for the truncated dynamics. 
Because the basis functions are related to fixed points of a renormalization group 
equation, Hamiltonians with fine scale degrees of freedom are equal to identical 
Hamiltonians with coarse scale degrees of freedom plus additional operators 
with additional fine scale degrees of freedom. Eliminating the fine scale degrees 
of freedom, rescaling and adjusting the parameters of the theory leads to a 
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renormalization group transformation. While the calculations are not trivial, 
implementation of the decoupling using similarity renormalization group method 
involves commutators of discrete canonical fields and algebraically computable 
coefficients. Methods for computing these coefficients were discussed in section 
five and implemented in the appendix for the parts of the free Hamiltonain that 
couple different scales. 

Another important property of the wavelet basis is that the scaling functions 
on any fixed scale are up, to an overall constant, a locally finite, compactly 
supported partition of unity. When inserted in the formal expressions for the 
Poincare generators, expressed as integrals of the energy momentum and angular 
momentum densities over a fixed time surface, the generator is decomposed into 
a sum of operators that act in different spatial volumes at a given time. This 
allows one to exactly test the Poincare commutation relations in finite volumes. 
These exact finite-volume components of the generators can be expressed in 
terms of the discrete wavelet fields by replacing each field that appears in the 
generator by its expansion in wavelet smeared fields. When these expansions 
are truncated, by eliminating small-scale degrees of freedom, the commutation 
relations are violated because the products of the small-scale degrees of freedom 
couple to the large-scale degrees of freedom. In the wavelet basis these correction 
can be identified and their relative importance can be calculated. 

A final important property of the wavelet basis is the ability to reduce all of 
the quadratures that arc needed in the theory to finite algebra. This requires re- 
placing conventional computational techniques that depend on functions looking 
smooth on small scales by new methods based on the renormalization group. 

We also demonstrate that SU(3) gauge invariance could be implemented ex- 
actly in wavelet truncated theories. While our discussion was limited to some 
illustrative topics, there are a number of other topics where the wavelet repre- 
sentation might have some advantages. These include the wavelet representation 
of the operator product expansion, wavelet representations of the Poincare Lie 
algebra in momentum space. 

This work was supported in part by the U.S. Department of Energy, under 
contract DE-FG02-86ER40286. The authors would also like to thank Profes- 
sor Robert Perry for useful discussions on the similarity renormalization group 
method and Andreas Schaefer for bringing our attention to some applications 
of wavelets in quantum field theory. 

References 

[1] I. Daubechies, Orthonormal bases of compactly supported wavelets, Comm. 
Pure Appl. Math. 41(1988)909. 

[2] I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Regional Conference 
Series in Applied Mathematics, 1992. 

[3] G. Kaiser, A Friendly Guide to Wavelets, Birkhauscr 1994. 



36 



[4] Howard L. Resnikoff and Raymond O. Wells, Jr.. Wavelet Analysis, Springer, 
1998. 



[5] O. Bratteli and P. Jorgensen, Wavelets through A Looking Glass - The World 
of the Spectrum, Birkhauser, 2002. 

[6 ] B. M. Kessler, G . L. Payne, W. N. Polyzou, Wavelet Notes, 
|arXiv:nucl-th/0305025|^2, 2003. 



[7] Christoph Best, Andreas Schaefer, "Variational description of statistical field 
theories using Daubechies' wavelets", arXiv: hep-lat/9402012 ( 1994 

[8] Christoph Best, Wavelet induced renormalization group for the Landau- 
Ginzburg model, Nucl. Phys. Proc. Suppl. 83 (2000) 848-850. 

[9] Ahmed E. Ismail, Gregory C. Rutledge, and George Stephanopoulos, "Mul- 
tiresolution analysis in statistical mechanics. I. Using wavelets to calculate 
thermodynamic properties" J. Chem. Phys. 118, 4414 (2003). 

[10] Ahmed E. Ismail, Gregory C. Rutledge, and George Stephanopoulos, "Mul- 
tircsolution analysis in statistical mechanics. II. 'The wavelet transform as 
a basis for Monte Carlo simulations on lattices" J. Chem. Phys. 118, 4424 
(2003). 

[11] Mikhail V. Altaisky, "Wavelet-Based Quantum Field Theory Symmetry, 
Integrability and Geometry: Methods and Applications", SIGMA 3 (2007), 
105. 

[12] P. Federbush, New formulation and regularization of gauged theories using 
a nonlinear wavelet expansion, Prog. Theor. Phys. 94 (1995) 1135-1146. 

[13] Guy Battle, Wavelets and Renormalization, World Scientific, Series in Ap- 
proximations and Decompositions, Volume 10, 1999. 

[14] S.Albeverio, M.V.Altaiskyb, "A remark on gauge invariance in wavelet- 
based quantum field theory" larXiv:0901."2806> 2. 2009. 

[15] E. Anderson, S.K. Bogner, R.J. Furnstahl, E.D. Jurgenson, R.J. Perry, A. 
Schwenk, Block Diagonalization using SRG Flow Equations, [arXii:0801. 1098, 
2008. 

[16] Okubo S., Prog. Theor. Phys., 1954, v.12, p.603. 

[17] Glockle W., Mueller L., Phys. Rev., 1981, v.C23, p.1183. 

[18] Kenneth G. Wilson and Wolfhart Zimmermann, Operator product expan- 
sions and composite field operators in the general framework of quantum field 
theory, Comm. Math. Phys. Volume 24, Number 2 (1972), 87-106. 



37 



10 Appendix - overlap integrals 



In this appendix we compute the overlap integrals D^ nn that appear in the 
free-field Hamiltonian. These are the terms responsible for the coupling of the 
degrees of freedom associated with different scales. 

The general coefficients are products of the matrices below in the x, y and 
z variables. There are nine combinations. 



D k mn = J a%(x)a%(x)dx = 2 2k J s'(x - m)s'(x - n) = 2 k D mn (233) 
D k ± = J a»{x)v%(x)dx = 2 2 (' +1 > H l ^7 k G nnl D m , nl (234) 

m'n' 

D&n = J wl' n (x)w l ^x)dx = 2 2 V+V (GH l - j ) mm ,G nn ,D m , n , (I > j). 

m'n' 

(235) 

The above expressions show that each of these integrals are linear combina- 
tions of the matrices 

D lm = ( a'^a'Ja^da. (236) 



We can use translational invariance to write (|236[) equation as 

Am = -Do.m-I = / s '{ x )s' m -i{x)dx. (237) 



Because of the support conditions on the scaling functions these vanish unless 
| n — m\ < 4. It is difficult to get an inhomogeneous equation for Di m because 
the obvious choice, using the partition of unity (|152j) satisfies 

J2 lD o,m-l=Q (238) 
i 

for any m. This is because what remains is the integral of the derivative of a 
continuous function with compact support. To get around this problem we use 
the partition of unity (|2ip to get the relation 



D , m -i = D lrn = ^2 D nim = ^2 J s n (x)s' l (x)s' m (x)dx. (239) 

n n 

This expresses the coefficients D mn as linear combinations of Di mn . Using the 
methods discussed in section five we find that coefficients D„; m = D 0l i- n ,m-n 
satisfy the Homogeneous equations 

D olm = 4V2 H o 

Hi,k+ n H mj+n D ok j (240) 

n,kj 

and 

J2 D oim - Doi m - Doml = (241) 
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D -4 


-5.357 x 


io- 


3 


Do-3 


-1.143 x 


io- 


1 


Do-2 


8.762 x 


io- 


1 


Do-i 


-3.390 






Doo 


5.268 






Doi 


-3.390 






D02 


8.762 x 


10" 


-1 


D 03 


-1.143 x 


10" 


-1 


Doi 


-5.357 x 


10" 


3 



Table 2: 



To-4 


-3.424658 x 10" 4 


To-3 


-1.461187 x 10" 2 


To-2 


1.452055 x 10" 1 


r -i 


-7.452055 x 10" 1 


Too 


-5.116622 x 10" 16 


r i 


7.452055 x 10" 1 


To2 


-1.452055 x 10" 1 


To3 


1.461187 x 10~ 2 


To4 


3.424658 x 10" 4 



Table 3: 



and the inhomogeneous equation 

^ IDnlm = / s n (x)s' m (x)dx — T nm — r , m -n- (242) 

I J 

An independent subset of these equation can be solved for the non-zero D i m s 
in terms of the To n . The coefficients To n satisfy the homogeneous equations 

To; = 2 Ho m Hi >n+m To n (243) 

m,n 

and 

J2 nT 0n = 1 (244) 

a 

where we have used the partition of unity (| 1521) . These equations can be solved 
for the non-zero To n and the solutions can be used to calculate Doi m from which 
one can derive D 0m using (|239l) . These can be used calculate (|233H235p and 
products of these quantities give the coefficients (|99M101[) for any combination 
of indicies. The results of the calculation of the nine non-vanishing Do m are 
given in Table 2. 

The results of the calculation of the non zero To n and Do m n are given in 
tables 3 and 4. 
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D 


_4_ 


4 


4.056756 x 


10 


-5 


-Do-1-2 


-6.544856 x 


10" 


1 


D012 


1.758631 x 10 


-1 


D 


-4- 


3 


1.620980 x 


10 


-4 


£>o-i-i 


2.323493 






A)13 


1.299066 x 10 


-2 


D 


-4- 


2 


-6.227505 x 


10 


-4 


-Do-io 


-2.071142 






-D014 


-8.819594 x 10 


-5 


Do 


-4- 


1 


9.394026 x 


10 


-4 


A>-ii 


3.284401 x 


10" 


-1 


-D02-2 


1.773552 x 10 


-3 


Do 


-40 




-5.193176 x 


10 


-4 


Dq-12 


3.304668 x 


10" 


2 


-D02-1 


3.304668 x 10 


-2 


Do 


-3- 


4 


1.620980 x 


10 


-4 


Dq-13 


-3.892382 x 


10" 


4 


D020 


-8.023148 x 10 


-2 


Do 


-3- 


3 


1.782152 x 


io- 


-2 


Aw -4 


-5.193176 x 


10" 


4 


D021 


1.758631 x 10 


-1 


Do 


-3- 


2 


-4.290543 x 


10 


-2 


-Doo-3 


-8.962795 x 


10" 


-3 


D022 


-1.066658 x 10 


-1 


Do 


-3- 


1 


4.009720 x 


io- 


-2 


-Doo-2 


5.753018 x 


10" 


1 


D023 


-2.400100 x 10 


-2 


Do 


-30 




-8.962795 x 


10 


-3 


A)o-i 


-2.071142 






Dq24 


2.149391 x 10 


-4 


Do 


-31 




-6.212589 x 


10 


-3 


-Dooo 


2.364229 






-D03-1 


-3.892382 x 10 


-4 


Do 


-2- 


4 


-6.227505 x 


10 


-4 


Adoi 


-7.734980 x 


10" 


-1 


-Do30 


-5.168265 x 10 


-3 


Do 


-2- 


3 


-4.290543 x 


10 


-2 


DQ02 


-8.023148 x 


10" 


2 


-D031 


1.299066 x 10 


-2 


Do 


-2- 


2 


2.549910 x 


10 


-1 


-D003 


-5.168265 x 


10" 


3 


£>032 


-2.400100 x 10 


-2 


Do 


-2- 


1 


-6.544856 x 


10 


-1 


DQ04 


-9.550006 x 


10" 


6 


-D033 


1.703740 x 10 


-2 


Do 


-20 




5.753018 x 


10 


-1 


An-3 


-6.212589 x 


10" 


-3 


-D034 


-4.695561 x 10 


-4 


Do 


-21 




-1.340525 x 


10 


-1 


-D01-2 


-1.340525 x 


10" 


-1 


D040 


-9.550006 x 10 


-6 


Do 


-22 




1.773552 x 


10 


-3 


An-i 


3.284401 x 


10" 


1 


Dq41 


-8.819594 x 10 


-5 


Do 


-1- 


4 


9.394026 x 


io- 


-4 


-D010 


-7.734980 x 


10- 


-1 


DoA2 


2.149391 x 10 


-4 


Do 


-1- 


3 


4.009720 x 


10 


-2 


Don 


3.965574 x 


10- 


1 


Dois, 
D044 


-4.695561 x 10 
3.523629 x 10 


-4 
-4 



Table 4: 
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